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When quantum back-reaction by fluctuations, correlations and higher moments of a state becomes 
strong, semiclassical quantum mechanics resembles a dynamical system with a high-dimensional 
phase space. Here, systematic computational methods to derive the dynamical equations including 
all quantum corrections to high order in the moments are introduced, together with a (deparameter- 
ized) quantum cosmological example to illustrate some implications. The results show, for instance, 
that the Gaussian form of an initial state is maintained only briefly, but that the evolving state 
settles down to a new characteristic shape afterwards. Remarkably, even in the regime of large 
high-order moments, we observe a strong convergence within all considered orders that supports the 
use of this effective approach. 
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I. INTRODUCTION 

Semiclassical approximations are of importance 
throughout physics in order to extract new effects or 
potentially observable phenomena in regimes in which 
quantum features are relevant but not dominant. All 
aspects of quantum gravity currently considered for po- 
tential tests fall in this class of situations. While deep 
quantum phases are crucial for several conceptual prob- 
lems, the current developments remain so diverse that 
reliable conclusions are difficult to draw. Semiclassical 
physics, on the other hand, is already useful for this 
class of theories in order to demonstrate their consis- 
tency with large-scale properties of the universe and to 
approach potential low-energy observations for instance 
via details of the cosmic microwave background. For this 
reason, the development of systematic semiclassical ap- 
proximations is one of the major requirements in current 
quantum-gravity research, especially for approaches such 
as loop quantum gravity whose underlying principles do 
not make use of a continuum geometry and perturbations 
around it. 

In loop quantum gravity and cosmology, the main tool 
available at a dynamical level is that of effective equa- 
tions which describe the evolution of expectation values 
in a dynamical state [l|. As in most quantum systems, 
especially interacting ones, quantum fluctuations, corre- 
lations and higher moments of the state then back-react 
on the evolution of expectation values, described in effec- 
tive equations by coupling terms between classical and 
quantum degrees of freedom. With an infinite number of 
independent moments of a state, there are infinitely many 
quantum degrees of freedom, making the general system 
of equations difficult to analyze. In semiclassical regimes, 
however, the infinite set of equations reduces to finite sets 
at any given order of the expansion by powers of h. In this 



way, a systematic expansion arises that goes well beyond 
other examples such as the WKB approximation. To first 
order in h and combined with a second-order adiabatic 
approximation, applied to anharmonic oscillators, equa- 
tions equivalent to those of the low-energy effective action 
are produced when expanding around the ground state 
0, Q. (The more widely applied WKB approximation 
does not reproduce the low-energy effective action [3(.) 
For other systems and different states, on the other hand, 
the canonical methods developed in |f| are more general 
and apply to quantum cosmology as well, even taking 
into account the totally constrained nature of these sys- 
tems [J]. One should thus expect that an application to 
quantum gravity allows comparisons with low-energy ef- 
fective quantum gravity but can also highlight new 
effects related to specific properties of quantum space- 
time structure as opposed to just graviton dynamics on 
a classical background. 

Compared with direct solution procedures of the par- 
tial differential equations that determine the dynamics 
of wave functions, the procedure used here has the ad- 
vantage of directly yielding expressions for expectation 
values and other, in principle observable quantities. The 
usual, but often numerically cost-intensive two-step pro- 
cedure of first finding expressions for the wave function 
and then performing integrations to obtain expectation 
values is thus significantly reduced. In this context, the 
methods used and described are useful not just for quan- 
tum cosmology, the realm of our examples, but may find 
wider applicability. In quantum cosmology, or more gen- 
erally constrained systems, other advantages include a 
possible extension to systems in which the dynamics is 
determined by a Hamiltonian constraint, not a Hamilto- 
nian function 0, 0, H| , and a simplification of the problem 
of time in semiclassical regimes ITo| . 

The semiclassical approximation allows one to con- 
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sider most cases of interest in current quantum-gravity 
research, but they do not allow one to tackle all of them. 
For this reason it is of interest to analyze the system 
of equations obtained as more and more of the mo- 
ments grow large and become relevant when a dynam- 
ical state turns more highly quantum. Since such large 
numbers of variables coupled to each other by long evo- 
lution equations are difficult to handle analytically (see 
the Appendix for deterrence), the main aim of this ar- 
ticle is to provide results of efficient computer-algebra 
codes designed to derive these equations automatically 
and quickly to high orders, which can then be fed into 
a numerical solver of the coupled differential equations. 
In this way, one can see how much high orders of the 
expansions matter and how higher moments can affect 
the evolution of expectation values as strong quantum 
regimes are approached. 

An interesting toy example for this is a spatially flat 
isotropic universe with a positive cosmological constant, 
filled with a homogeneous free and massless scalar field. 
The free scalar can play the role of a global internal time, 
allowing one to reformulate the constrained system as or- 
dinary Hamiltonian evolution. In the presence of a posi- 
tive cosmological constant, the value of the scalar field is 
bounded from above as a function of time; thus, when us- 
ing the scalar as internal time, infinite volume is reached 
at a finite time 4>div As this point is approached, volume 
fluctuations and higher moments diverge even though the 
regime is supposed to be one of low-curvature nature (see 
Figure [T|) . The non-semiclassical appearance is an arti- 
fact of the choice of time, highlighting the toy-nature of 
the model. Physical conclusions derived from the model 
for a late-time universe with a positive cosmological con- 
stant can thus not be reliable, as interesting as they might 
be. But the model is ideal for a mathematical analysis 
of the possible roles of large moments. In the present 
article, we will use this example mainly to illustrate the 
usefulness of the numerical codes; further analysis re- 
quires more sophisticated methods to handle the large 
parameter space of all moments. 



II. GENERAL FORMALISM 

We begin with a brief review of the procedure of effec- 
tive equations, followed by a correction of an important 
formula of [1] and its accompanying proof. 

A. Quantum state space 

Classically, the phase space of a system with one degree 
of freedom is completely described by the generalized po- 
sition V (related to the volume in cosmological models) 
and its conjugate momentum P (related to the Hubble 
parameter). During the quantization process each clas- 
sical phase-space function is promoted to an operator, 
and observable information is extracted via expectation 
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FIG. 1: In this plot we show the evolution of the volume 
expectation value V for a cosmological constant A = 9 x 10 7 
and initial data to be presented in Sec. lIVEl together with the 
spread around this trajectory V ± AV. This plot illustrates 
the growing dispersion of a state evolving to large volume. 
The range for <j> shown, relative to the value 4>&iv at which the 
classical volume diverges, is the maximal one obtained with 
numerical stability to tenth order of moments of the evolving 
state. As the volume grows, the area swept out by the wave 
packet more and more aligns itself with the vertical axis, such 
that volume fluctuations (as well as some other moments) 
computed at a fixed value of <f) diverge. This divergence would 
happen even if the state remained sharply peaked transversal 
to the direction in which it moves, and is thus not a failure 
of semiclassicality in this system but rather an artifact of the 
parametrization of evolution. Physical implications are thus 
difficult to draw, but an interesting test case for back-reaction 
issues is obtained. 



values which can be seen as functionals on the algebra of 
operators that characterize states. Since expectation val- 
ues of products of operators in general differ from prod- 
ucts of expectation values, the following infinitely many 
moments are needed to describe the system completely: 

G a ' b := ((P - Pf (V - V) b ) Wcyh (1) 

where the subscript Weyl stands for totally symmetric 
ordering. Moments cannot be chosen arbitrarily but are 
subject to additional conditions such as the uncertainty 
relation 

G 2 '°G ' 2 - (G 1,1 ) 2 > ^ . (2) 

Moreover, a set of moments satisfying all those condi- 
tions in general corresponds to a mixed rather than pure 
state. By considering moments to describe states, mix- 
ing (which is important for generic descriptions of ho- 
mogeneous quantum cosmology seen as an averaged de- 
scription of microscopic physics, in particular inhomo- 
geneities) is automatically included. 

The use of moments allows a very general definition 
of semiclassical states as those satisfying the hierarchy 
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G a ' b = 0(ft( a+b '/ 2 ), realized for Gaussian states (see be- 
low) but also for a much wider class. This character- 
ization of semiclassical states (pure or mixed) has two 
major advantages: (i) It is a general characterization, as 
unbiased as possible. (In more restrictive classes such 
as Gaussians the specific form of the wave function may 
actually matter for physical effects.) (ii) The infinitely 
many G a,b are decomposed in finite sets to order H n . This 
feature allows systematic approximations of semiclassical 
dynamics. 

The moments are pure quantum degrees of freedom; 
they would all vanish in a classical limit. In a quan- 
tized interacting system they couple to the expectation 
values P := (P) and V := (V). The state with respect 
to which these moments are taken is specified via values 
for all the quantum variables, which are all dynamical 
and satisfy equations of motion. In order to consider the 
evolution of expectation values, the evolution of all mo- 
ments must be followed simultaneously. Assuming the 
moments to be constant or to take specific values such 
as, for instance, those of a Gaussian state at all times 
cannot provide the correct dynamics. The dynamics is 
described by a quantum Hamiltonian, the expectation 
value of the Hamiltonian operator, which is state depen- 
dent and thus a function of expectation values and all 
the moments. By Taylor expansion, one can write the 
quantum Hamiltonian as 



is still a corresponding wave function with those moments 
to within the order of approximation considered. It may 
be difficult to construct such a wave function explicitly, 
but this is not required because the information relevant 
for observations is already contained in the expectation 
values and moments. (When moments diverge, as in the 
example discussed below, the state moves out of the do- 
main of definition of the operators considered. By our 
moment expansions we will only consider the approach 
to such a point of divergence, not the point itself.) 

In order to obtain all Hamiltonian equations of mo- 
tion in the context of this effective approach, we will 
need to compute the Poisson brackets between any two 
moments. If one defines the Poisson bracket for expecta- 
tion values of arbitrary operators X and Y by the rela- 
tion {(X), (Y)} = -ih~ l {[X,Y]) (extended to arbitrary 
functions of expectation values using linearity and the 
Leibniz rule), Hamiltonian equations of motion generated 
by Hq = (H) are equivalent to the Schrddinger flow of 
states or the Heisenberg flow of operators. By applying 
this relationship to all powers of the basic operators, the 
Poisson brackets between moments is obtained. Since 
[V, P] — ih, in the particular case of the basic operators 
themselves, the Poisson bracket reduces to the classical 
one, {V,P} = I. Moreover, one can easily show that 
moments have vanishing Poisson brackets with the basic 
expectation values: 
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with the classical Hamiltonian H. The quantum vari- 
ables it contains provide crucial corrections to the clas- 
sical dynamics in a state-dependent way. We note that 
expansions of this form (as also the Feynman or WKB ex- 
pansions) are expected to be asymptotic, not necessarily 
convergent. In the model discussed below, we will nev- 
ertheless observe strong convergence properties within a 
large range of orders. 

The evolution generated by the quantum Hamilto- 
nian Hq , governed by the equations of motion d/ /d(f> = 
{/,Hq} for phase-space functions /, is described with 
respect to a time parameter <p. This could be an abso- 
lute time parameter or, as in the example used below, 
an internal time variable chosen for a deparameterizable 
constrained system. (For the effective treatment of non- 
deparameterizable constrained systems, lacking a global 
internal time, see [13].) Effective quantum evolution then 
provides an approximation scheme for time-dependent 
expectation values of operators, starting with some ini- 
tial state specified by the initial expectation values and 
moments. It is straightforward to find moments corre- 
sponding to a wave-function representation of the initial 
state in a Hilbert space, such as the Schrddinger one. Mo- 
ments evolved by the quantum Hamiltonian then corre- 
spond to the dynamical state, evolved in the Schrddinger 
or Heisenberg picture. We are thus guaranteed that at all 
times, as long as the evolved moments remain finite, there 



{G a ' b ,V} = Q = {G a < b ,P} 



(4) 



The general case of Poisson brackets between moments 
is more involved. We have the general formula 



{G a ' b ,G c > d } = adG a - 1 ' b G c ' d - 1 - bcG a ' b - 1 G c ~ l '' 

n-l 



E( ih 
hr 



l^n s~ia-\-c—n,b-\-d—n 
IX abcd U ) 



(5) 



where the sum over n runs over all odd numbers from 1 to 
Min(a + c, b + d, a + b, c + d), which makes the coefficients 
real in spite of the presence of the imaginary unit. We 
have also defined 
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where the C-coefficients come from the Weyl ordering 
(see also Sec. IIII[) and are defined as 



C\ 



i\k\ 
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d\ 



(n - d)\(k - d)\(2d)U 2 d 



k 



(7) 



Note that these coefficients have the permutation prop- 
erty K™ bcd = (— l) n K™ dab . Since the sum in §5$ is only 
over odd n, this makes the antisymmetry of the Pois- 
son brackets transparent. (In these formulas, G 0,0 = 1, 
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£0,1 = G i,o = Q) and G a,b _ o for any a or 6 < are 
understood.) 

A general formula analogous to (J5J) was obtained in [l| , 
but since it was not applied to moments of high orders, 
a typo remained undiscovered. This typo is corrected in 
©, and since the formula plays an important role for the 
results of this article, we provide a detailed proof in the 
following subsection. 



B. Poisson structure of moments 

As in we consider the generating function, 
D(a) = ( e «p(*-P)+av(fr-V)) 
1 (3 
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G j - k ' k ,(12) 



ml(j — k — m)\n\(k — n)\ 
whereas, for the second one, the expansion is given by 
(ap[3 v - a v j3p)D{a)D{j3) 

oo j oo m /~yj—k,kf~im—n+n 

^ — > ^ — "V ^ — > ^ — > Kjr Lx 

2— 1 2— 1 2— 1 2— 1 k\(j - k)lnl(m - n)\ 

j=0 k=0 m=0 n=0 w ' v y 

x (a J p - fc+1 at>/??^"/3^ +1 - a^ fe a^ +1 /3™"™ +1 /3^).(13) 

Since ([TTt must be satisfied for all values of a and j3 we 
find that the coefficients in front of each of the powers of 
various a's and /3's must be equal. By direct comparison 
we then obtain the Poisson bracket between our quantum 
variables. Making the convenient substitutions a — w—x, 
b = x, c = y — z, and d = z, and combining the factorials 
into binomial coefficients the result reads, 



Assuming that the wave function is analytic in the mean 
values, D(a) and the Poisson bracket between two of 
these functions will also be analytic, which allows us to 
take the Taylor expansion to all orders in a, 



oo w oo y 



w P - x a* v ri P - z p* v 



w=Q x=0 y=0 z=0 
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x\{w — x)\z\(y — z)\ 
(9) 



On the other hand, we use the Baker-Campbell- 
Hausdorff formula to get the commutator 
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which gives us the Poisson bracket 



{D(a),D((3)} = -sin \hotv n r h ) : > M 



(a v (3p-a P (3v)D(a)D((3). 



(11) 



We then use a Taylor series and the binomial theorem to 
expand the right-hand side of the last equation. For the 
first term we get the expression, 

| sin (-(a v j3 P - app v yj D(a + (3) 
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(14) 



Finally we make the replacement n = 2r + 1 , noting that 
in the summation over r the only non-zero terms occur 
when 2r + 1 < Min(a + c,b + d, a + b,c + d), which gives 
us ©. 



III. ALGEBRAIC IMPLEMENTATION 

We have written two different codes in Mathematica 
in order to obtain the Poisson brackets between any 
two generic moments. The first code works iteratively, 
whereas the second one makes use of the general formula 
(O. The results agree for all orders checked (up to 14th) 
but the second code is much more efficient. This can be 
clearly seen in Fig. [2] where we have plotted the timing 
of both codes to compute several Poisson brackets. The 
timing of the recursive code increases exponentially with 
the order of the considered Poisson bracket, whereas the 
general code allows us to obtain high-order results in a 
very short time. This fact shows the importance of the 
general formula (J5]) from a practical point of view. 

In the following we briefly explain how the recursive 
code works. Making use of the binomial identity, the 
expression for the moments can be written as, 



Ga - b = EEt- 1 ) 



+b-k-n 



k=0n=0 

y^pa—k-yb—n ipk y-rn 



Weyl- 



(15) 
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FIG. 2: In this figure we compare the timing of computing 
the Poisson brackets between any two moments with the code 
working iteratively as well as the one that makes use of the 
general formula ([5]). The time, shown in a logarithmic scale 
(to base 10) , has been measured in seconds. The order plotted 
in the graphics corresponding to a bracket {G a ' b , G c ' d } has 
been defined as the sum of all the indices: a + b + c + d. Note 
the exponential increase of the timing for the recursive code, 
whereas the computations with the general code are kept in 
the order of 10 -4 seconds up to the considered order. There 
are N = (n + 1) (n + 2)/2 — 3 different moments with an order 
less or equal to n, where the subtraction of 3 comes from the 
fact that the order and 1 moments are trivial; they can be 
combined into N(N + l)/2 different Poisson brackets. In the 
case of the recursive (general) code, the 33 (228) moments 
up to order seven (20) have been considered, and in the plot 
the timing for their corresponding 561 (26106) independent 
Poisson brackets are drawn. 



On the other hand, we also have found that the Weyl 
ordering can be given as a linear combination of expec- 
tation values with a predefined order, for any operators 
A and B, 

Min(fc.n) 

(i fe S")wcyi= £ [B,A] d Ci n (A k - d B n - d ), (16) 

d=0 

where the coefficients C d n have been defined in ([7]). 
Hence, for our case of interest, we can write either 

Min(fc,n) 

(p fc o Wcyl = m d cL(p k ~ d v n - d ), (17) 

d=0 

or. 

Min(fe,n) 

(P k v n ) Wcyl = ]T (-ih) d c d kn (v n - d P k - d ). (is) 

d=0 

In order to obtain the Poisson brackets between different 
moments iteratively we make use of the definition of the 
moments in the form, 

Ga * = EE E {-m d (~iy +b - k - n ( a k )( n ) 



xC d n P a - k V b ~ n {V n - d P k - d ). (19) 

For computing {G a ' b , G c ' d } we then use the properties of 
the Poisson brackets as well as of the commutator until all 
operator terms reduce to the simple commutator between 
the isolated operators V and P. 



IV. A MASSLESS SCALAR FIELD WITH 
COSMOLOGICAL CONSTANT 

As an application of high-order quantum back-reaction 
we now introduce a cosmological model with a regime of 
large moments: an isotropic spatially flat universe filled 
with a free, massless scalar <p (with momentum p^) and 
a positive cosmological constant A > 0. From the Fried- 
mann equation 

we obtain, ignoring constant factors for the sake of com- 
pactness, = a 2 y/a' 2 — Aa 2 where the prime is the 
derivative with respect to cosmological time. The square 
root may take both signs, corresponding to expanding 
and contracting cosmological solutions. We will choose 
the positive sign in order to analyze expanding universe 
models approaching large volume, as they remain in a 
low-curvature regime. The opposite sign, if desired, could 
be easily included since it would just invert the sign of the 
constant of motion p$ [which will be the deparameterized 
classical Hamiltonian of our system (f2"2j)] as well as the 
effective Hamiltonian (|2U| below. Therefore, from the 
expanding solutions we show it will be straightforward 
to get their contracting counterparts just by performing 
a time reversal. 
We rewrite p^ as 

p = (1 - x)V^P 2 - A[(l - x )v]0-+2x)/o-x) (2i) 

in terms of canonical gravitational variables V = (1 — 
x) _1 a 2 ^ 2:c and P = —a 2x a'. The parameter x charac- 
terizes different cases of lattice refinement of an underly- 
ing discrete state; see [Tl[ [l2j . The quantum dynamics 
(Wheeler-DeWitt type) that will be implemented in this 
paper is insensitive to the change of this discreteness pa- 
rameter. However, some values lead to advantages in the 
solution procedure, as we will make use of below. More- 
over, it is useful to keep the parameter in the definition 
of canonical variables in order to facilitate possible ap- 
plications to loop quantum cosmology in the near future, 
where this parameter will become relevant. 

The system is simplest to deal with for x = —1/2, in 
which case V is proportional to the volume, because p^,, 
which can be thought of as the Hamiltonian for evolution 
in internal time </>, is then linear in V. For a negative cos- 
mological constant, the effective evolution coming from a 
quantization of (f2"Tj) has been solved and analyzed in [lj[ , 
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some of whose formulas apply here as well. (For numer- 
ical wave-function evolution of this model with A < 0, 
see [HI].) For A < 0, however, the moments do not grow 
large at large volume, and the volume is bounded by rec- 
ollapse. These two features are not present for A > 0, 
making an analysis of high-order moments necessary, the 
topic of this article. 

From now on we use x = —1/2, which implies 



= H = -F\/P 2 -A 
2 



(22) 



as the Hamiltonian of the deparameterized system. Since 
the dependence of the classical Hamiltonian on V is lin- 
ear, we will only obtain terms of the form G a '° and G a>1 
in the quantum Hamiltonian and we have the closed ex- 
pression, 



Equations of motion for the quantum variables can also 
be derived from Poisson relations between them 



3 ^™ A-"/ 2 

n=2 



VT n {P/VA){G r ' s ,G n > } 



+WAT„_ 1 (P/v / A) {G r ' s , G™ -1,1 } 



.(28) 



In fact, this highly coupled system of infinitely many vari- 
ables is extremely complicated to solve, and we will try 
to analyze it by means of truncations to finitely many 
variables which include the expectation values and only 
low order quantum variables. 



q q x ' A-n/2 

h q = |^v^a + |vaeV 



n=2 

,VAT n _ 1 (P/VX)G n - 1 ' 1 



VT n (P/VA)G n '° 
(23) 



where we have defined 



A. Vanishing A 

For A = 0, we obtain the harmonic model of [l5[ for 
which all equations decouple: 



Explicitly, the moments appearing here are 



(24) 



G h 



((p-(p)T), 



1 n— 1 

G n-i,i := _j2((P-(P)y(v-(v))(P-(P)r- 1 - 



Hamiltonian (|23l) will describe the semiclassical behavior 
of the Wheeler-DeWitt quantization of this system. 

Equations of motion for V and P are derived as in any 
Hamiltonian system, noting that the expectation values 
Poisson commute with any quantum variable. This gives 



V 



P 



2 V VP^A ' 2 ~ „ 



3g A 1 



VT n+1 (P/VA)G n >° 



+ nVAT n (P/VA)G 



V 
P 



G ab 
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(29) 
(30) 
(31) 



There is no quantum back-reaction, and all equations can 
easily be solved: 



v(4>) 
p(4>) 

G ab (0) 



Vbexp 
P exp 



Gg h exp 



(25) 



-A4>-<M 
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(32) 
(33) 
(34) 



P 



Q o 00 A-n/2 

VP 2 -A- -x/A^T — — T n (P/VA)G"'°, (26) 



n=2 



where dot stands for derivatives with respect to </> and 
we have made use of the fact that 



[v,T n (P/VX)} = -±=T n+1 (P/VA) 



(27) 



B. Classical equations 



To classical order, i.e., assuming all moments to vanish, 
we obtain analytical solutions 



P 



lassical 



(0) = P cosh(3(0 - 0o)/2) - ^P 2 -Asinh(3(<£ - o )/2) , 



(35) 



7 



Vr 



classical 



(</*) = V 



A 



v/P 2 -Acosh(3(^ - 0o)/2) - P> sinh(3(<^ - o )/2) ' 



(36) 



r 



where Vq and Po are the value of the expectation val- 
ues at the initial time </>o. At the time when <p^ v := 
0o+2/3tanh~ 1 ( v /P o 2 - A/P ) the volume Kl aS sical(^) di- 
verges, and v/Pclassical^div) 2 - A = 0. 



C. Second order 

To understand whether quantum effects can be impor- 
tant in this large-volume regime of small curvature (if 
A is small), we will analytically solve (under certain ap- 
proximations) the equations of motion for second order 
quantum variables, i.e. fluctuations G 2,0 , G 0,2 and the 
covariance G ' . To this order, ignoring higher moments, 
we have a quantum Hamiltonian 



V 



H * = 2^^-lV-A)3/ 2 



G 



2.0 



P 



--G 



i,i 



2 VP 2 - A 
Using the Poisson bracket relations 

{G°< 2 , G 1 * 1 } = 2G°< 2 , {G°< 2 , G 2 <°} = 4G 1 * 1 
{G^.G 2 - } = 2G 2 ' , 

the moments satisfy the equations of motion, 
P 



(37) 



(38) 



G 2 < u = 



G 1 ' 1 



VP 2 - A 

3 V 

"2 (P 2 



=G 2 



A) 3 /' 



-G 



2.0 



G ' = -3A 



r 



(P 2 - A) 3 / 2 



G 



VP 2 - A 



(39) 
(40) 
G°< 2 .(41) 



i. Expectation values as background solutions 

As a first approximation, we can use the classical solu- 
tions for V and P as some kind of background on which 
the quantum variables evolve. From the resulting solu- 
tions we then determine whether quantum back-reaction 
effects by quantum variables on expectation values are 
significant. 

The equations for fluctuations to this level of trunca- 
tion are not strongly coupled, which allows us to solve 
for all of them by integrations. For G 2,0 , we obtain 

G 2 '°(0) = G 2 ^(0o)fcosh(3(^~^o)/2) 



P, 



sinh(3(0^ o )/2) , (42) 



which then allows us to solve for G 1,1 and, in turn, G ' 2 . 
We do not present these solutions here (which can be 
found from [13] by changing the sign of A), but note that 
fill) can be written as 



G° 



V 2 

classical 



= -3A- 



G 1 



Vr 



classical 



(P 2 - A) 3 / 2 



(43) 



(Within the second-order approximation, we could re- 
place V^iassicai with V in this equation up to terms of the 
form G°- 2 G 2 ' and G°' 2 G 1 ' 1 .) This shows that the evo- 
lution of volume fluctuations relative to volume is deter- 
mined by the covariance G ' . For an unsqueezed state, 
we would have G 1,1 =0 and squared volume fluctuations 
would be proportional to the total volume. However, 
even if we start with an unsqueezed initial state, squeez- 
ing would develop over time since G 1 ' 1 is non-zero. In 
this way, the precise state has to be analyzed to under- 
stand the long-term evolution of fluctuations. 



2. Coupled equations 

The solution (|42|) for G 2: °, compared with (|36|) . shows 
that curvature fluctuations to second order vanish where 
V diverges. Thus, also G ' 2 has to diverge there in or- 
der to respect the uncertainty relation. The approach 
to zero is the same for v G 2,0 , 1/V and VP 2 — A, which 
we can use to check whether quantum back-reaction is 
important. Writing equations of motion for expectation 
values to the order where coupling terms to fluctuations 
appear, we obtain 



P 
V 



G 



2.0 



VP 



2 VP 2 - A 



9 VPG 2 ' 

4 (P 2 - A) 5 /2 



-A 



G 1 



2 (P 2 - A) 3 / 2 ' 



(44) 



(45) 



which now forms a coupled set of five differential equa- 
tions if we include G° :2 . When V diverges, coupling 
terms due to quantum back-reaction cannot be ignored. 
In P, for instance, the coupling term diverges, and in 
fact dominates the classical term, since G 2,0 vanishes 
like P 2 — A, and also V has diverging coupling terms. 
Even the truncation used here, which does include some 
quantum corrections, is not consistent at this point be- 
cause one has to expect that further correction terms by 
higher moments become large, too. Nevertheless, in the 
approach to the divergence of V, correction terms should 
become relevant one by one, depending on their moment 
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order. Thus, it is of interest to analyze the coupled sys- 
tem written here to second order in the moments before 
proceeding to higher orders. 

We first focus on the system given by the equations 
(|44f for P and ([39]) for G 2,0 since these two equations 
are coupled only with each other. If we divide P by G 2 ' , 
we obtain the differential equation 



IP - A 1 



A 



dP 

dG 2,o - 2 PG 2 <° ~~ 4 P(P 2 - A) ' 

for P(G 2 ' ) or, in a simpler form, 

d(P 2 -A)_ p2 _ A _l_A_ G2 , . 



dlogG 2 ' 
This equation of the form 



2P 2 - A 



(46) 



(47) 



f , (x) = f(x)-lA- ea 



2 f{x) 



can be solved by 



f(x) = \/ce 2x + Ae x 

with an integration constant c. In this way, we derive the 
exact second-order relations 



P 



• AG-'". (48) 
v/P 2 - A = </c(G 2 > ) 2 + AG 2 ' . (49) 

We can use these relations to solve the equations of 
motion for P(</>) and G 2,0 ((f>), but the resulting differen- 
tial equations are complicated. For G 2 ' , for instance, 
the differential equation is 



G 2 ' 



P 



VP 2 - A 



G 



2,0 



-3,/(G 2 >°) 2 + 



AG 2 ' 



+ A/G 2 ' 
(50) 

Fortunately, as we will see in the next subsection, al- 
ready the relation P(G 2 ' ) has interesting implications 
regarding the quantum back-reaction problem. We will 
continue using (|50|) in the next subsection, but first note 
that further analysis, at least numerically, of the trun- 
cated second order system could be facilitated by an 
explicit decoupling of the equations. We have already 
provided direct equations for P and G 2 ' which we can 
assume to be solved by integrations. Further equations 
are then still coupled amongst each other, but different 
combinations are decoupled. First, we combine V and 
G 1 ' 1 to 



d_ /G 1 - 1 



V 



3 . G 2 ' 
-A 



P 



G 



i,i 



2 (P 2 -A) 3 / 2 2 VP 2 - A V 
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A- 



P 



4 (P 2 -A) 5 /2^ V 
3 A 1 ^ un2 



:G 2 <° 
G 



G 



i,i 



2 (P 2 - A) 3 / 2 V V 



(51) 



which is already decoupled from the rest. Once G 1,1 /V 
has been solved for, ([4"5"j) presents an uncoupled differen- 
tial equation for V, and we can finally integrate (|41|) for 
G 0,2 . 



3. Strength of quantum back-reaction 

The term AG 2 ' in (|4T)|) arises through the quantum 
back-reaction in the evolution equation for P ([44]) . If 
this term is ignored, we reproduce the previous approxi- 
mated solution G 2 ' = G 2 '°(</>o)(P 2 - A)/(Po- A), which 
follows from ([35]) and ([42]). Hence, comparison shows 
that the integration constant is proportional to c ~ 
(Pq — A) 2 /[G 2: °((/>o)] 2 if initial values are imposed where 
back-reaction is weak. Starting with a semiclassical ini- 
tial state, for some time we will have G 2,0 (cf>) ~ G 2 '°((f>o) 
and the term AG 2 ' can indeed be ignored compared to 
the classical parameter (Pq — A) 2 in (T4T)]) . (The value 
of Pg — A cannot vanish under the stated assumptions, 
for we already know that back-reaction terms would then 
be large at 4>q.) But as we approach the divergence of 
V, we have seen that G 2 ' according to the lowest or- 
der solution (|4"2"|) would tend to zero. In this regime, the 
new term A/G 2 ' in (|50[) will eventually dominate over 
the constant c. This is another demonstration of the fact 
that quantum back-reaction is important at this place. 

For c = 0, i.e. in the strong quantum back-reaction 
regime where A/G 2 ' dominates, we can find a solution 
for (E 



The combination h{4>) := ^/G 2 '°((/))/A satis- 
fies a simple differential equation solved by G 2,0 (</>) = 
Asinh 4 (3(0— (f>o)/4). The higher power compared to the 
solution (|42|) in the absence of back-reaction again indi- 
cates that quantum back-reaction becomes important at 
large volume. 

Turning now to V we can confirm that the quantum 
back-reaction term does change its behavior of diver- 
gence. In (|45p. the dominant back-reaction term involves 
G 2 ' rather than G ' . Ignoring quantum back-reaction 
in the solution for G 2,0 , we would have 



V 
V 



P 



AP 



2VP J3 A 4^(P 2 -A)3/2 



The second term is subdominant. Quantum back- 
reaction in the solution for G 2,0 , on the other hand, im- 
plies a stronger divergence as can be seen by using rela- 
tion (14"81 for c = 0, since we now have 



V 
V 



15 



P 



4 VP 2 - A 



Here, the quantum back-reaction term for V has the 
same form as the classical term, which have thus be com- 
bined. The pre-factor of the right-hand side becomes 
larger than classically once quantum back-reaction is in- 
cluded, strengthening the divergence of V. 

One may conjecture that, with all quantum corrections 
from the complete series of moments, the divergence of 
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V would disappear. This is suggested by the fact that 
a self-adjoint Hamiltonian would provide unitary evo- 
lution in which the wave function remains well-defined 
throughout the point where the classical V diverges, and 
even beyond that point. In the model considered, the 
Hamiltonian for 0-evolution is not essentially self-adjoint 
but has several inequivalent self-adjoint extensions [l6| . 
While subtle issues of self-adjoint extensions of operators 
usually do not play a large role in effective equations, pro- 
vided that consistent truncations are possible, the point 
where V diverges classically requires all moments to be 
taken into account. With infinitely many independent 
variables, coupling terms have to be arranged in a pre- 
cise manner, which can reflect issues of self-adjoint exten- 
sions. When all moments are taken into account in a way 
corresponding to unitary evolution of the corresponding 
wave function, well-defined equations for expectation val- 
ues and the moments may result. This supports the con- 
jecture that the classical divergence might be removed 
by an interplay of all moments, but a verification in this 
highly coupled system is complicated. (It could, in fact, 
happen that wave functions evolve out of the domain of 
definition of the unbounded operator V , and then di- 
rectly re-enter. Thus, it is not obvious that equations 
for expectation values (V) must be well-defined at all 
times. Another example which indicates caution against 
applying arguments of self-adjointness for the behavior 
of expectation values is the upside-down harmonic oscil- 
lator. Its Hamiltonian is not essentially self-adjoint and 
allows different self-adjoint extensions, depending on the 
boundary conditions at infinity where wave packets are 
reflected. However, the Hamiltonian is also quadratic, 
and no quantum back-reaction results. Effective equa- 
tions are automatically truncated without implementing 
any approximation; the evolution of moments to finite or- 
der just does not depend on what self-adjoint extension 
is used.) 

We will probe the divergence of V further by the nu- 
merics of high-order moments. For now, conclusions one 
can draw already from the analysis presented so far are: 

• Quantum back-reaction is essential at large volume 
in the presence of a positive cosmological constant, 
even though one would expect it to be a semiclas- 
sical regime of small curvature if A is small. 

• While a wave packet in a ^-representation may 
not show large deviations from the classical tra- 
jectory, this is only because volume fluctuations 
grow and are larger than deviations from the clas- 
sical trajectory. Nevertheless, as volume fluctua- 
tions diverge, deviations from the classical trajec- 
tory can be large. Equations including quantum 
back-reaction clearly show that correction terms 
cannot be ignored. Quantum back-reaction seems 
to strengthen the divergence of V rather than trig- 
gering a recollapse. 

• Corrections are more visible for curvature P since 
its fluctuations approach zero while also here quan- 



tum back-reaction terms in the equations of motion 
are relevant. Thus, deviations from the classical 
trajectory build up and, unlike for V, are not cov- 
ered by a broad wave function if a P-rcprcsentation 
is used. 

• All state parameters matter for the spreading be- 
havior and for the exact size of quantum correc- 
tions. A state may be assumed to be unsqueezed 
initially, setting G 1,1 (0o) = 0, but correlations will 
build up over time. 

D. Third order 

At third or higher orders, the equations of motion not 
only become longer; there is also a new feature in the 
Poisson relations which we illustrate here. Some exam- 
ples of third-order moments are 

G 3,0 = (p3) _ 3p G 2,o _ p3 ) 

G 2 ' 1 = -{VP 2 + PVP + P 2 V) -2P G M 

- V G 2 ' - P 2 V, 

G 1 ' 2 = -{PV 2 + VPV + V 2 P)~2VG 1 ' 1 

- P G°< 2 - V 2 P. 

The Poisson brackets of second-order with third-order 
moments are of third-order form, for instance: 

{G°< 2 , G 3 < } = 6G 2 * 1 , {G - 2 , G 2 * 1 } = 4G 1 * 2 , 
{G M , G 3 ' } = 3G 3,0 , {G X ' X ,G 2 ' X } = G 2,1 , 
{G 2 <°,G 2 < X } = -2G 3 <° , {G 2 <°,G 3 < } =0, 

where only the ones needed for evolution in our cosmo- 
logical system were written. Two third-order moments 
in a Poisson bracket, on the other hand, show a different 
behavior: They produce a term of higher order (fourth) 
as well as products of second-order moments, which are 
to be considered as being of fourth order in the hierarchy. 
This behavior of the Poisson structure has consequences 
for the truncation of equations of motion. 

Unlike for second order, (|38| . there is a difference be- 
tween expanding the Hamiltonian to third order and then 
deriving equations of motion, and expanding the equa- 
tions of motion (obtained from an expanded Hamilto- 
nian or the whole series in h) to third order. In particu- 
lar, moments of third order in the Hamiltonian produce 
moments of fourth order in the equations of motion for 
third order moments. This would result in a complete set 
if also equations of motion for fourth order moments are 
included. This problem does not arise if one first com- 
putes equations of motion and expands those to third (or 
any) order which is anyway more suitable because it is 
the equations of motion which we are primarily interested 
in. We illustrate this feature by the following evolution 
equations, computed from the third-order Hamiltonian: 
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p - ^^^ + 4 A (P^-4V^ G3 ' ' ( 52 ) 

V = 3 VP 9 a2 VPG 2 ' 3 2 GM 3 2 y(4P 2 + A) 3.o , 9 2 P 2>1 

2 VP 2 - A 4 (P 2 -A) 5 /2 2 (P 2 -A) 3 / 2 4 (P 2 -A) 7 /2^ ^4 (p2_A)V2° 

G 2,0 = -3 - f G 2 -° + H — ry-G 3,0 , (54) 

VP 2 - A 4 (p2_ A )3/2 

rU _ 3 AF 2;0 9 APT/ 3;0 3 A 2;1 

2(P2-A) 3 / 2 4(P 2 -A) 5 / 2 4(P2-A) 3 /2 ' 1 j 

^0.2 = o Ay u P . 2 9 AFP 3|0 3 A 2<1 

(P2-A)3/2^ VP 2 -A 4(P 2 -A)5/2 LT 4(p2_ A )3/2^ ' 

(56) 

I 

If we truncate the Hamiltonian (and hence the equa- shown in the following equations to third order in mo- 
tions of motion) to a given order, the equations of motion ments, 
for quantum variables involve higher order moments, as 



G3 '° - -f vw=t G "° + 1 ((G2 ' 0)2 - G4 ' 0) ' (57) 

~ ~2 A (P2-A)3/2 G + 4(P2-A)5/2 } ~ G j + 2VP 2 -A ' (58) 

G1,2 = _3A (P2-A)3/2 G2,1 - 2(P2-A)5/2 ( Gl,lG2 '° ~ G3,1 ) + 2VP2-A G2>1 ^ 

-i (p2 _ A A) 3/2 (-4(GM)2 +G 2,0 G 0.2 +3G 2, 2) 

AO, 3 _ 9 V x 2 27 AFP ^,2^2.0 r*2,2\ I ^ P q 3 

~ "2(P2-A)3/2 G -T(P2-A)V2 1 G - G ^ + 2VP2-A 

4 (P2 A A )3/2 ( GMG °' 2 - G " 3 ) ■ ( 60 ) 



For consistency of the approximation one should elimi- 
nate in this third-order approximation all terms of quar- 
tic order, i.e. fourth-order moments, producing a closed 
set of equations. (One may eliminate terms quadratic in 
second-order moments as well, but they are anyway sub- 
dominant when the hierarchy is satisfied.) Truncating 
only the Hamiltonian to third order is thus not enough 
because the general Poisson bracket between two third- 
order moments produces terms of fourth order. (In gen- 
eral, the Poisson bracket of two moments of order n and 
m, respectively, produces terms of order n+m—2 [see Eq. 
(J3J)]. Only for second order does this result in a closed 
Poisson algebra.) Consistently truncated equations of 
motion are produced if one truncates the Hamiltonian as 
well as the Poisson brackets to the order considered. An 
example to order five can be found in the Appendix. 



E. High-order corrections 

Because of the complicated structure of the equations 
at higher orders, analyzing them by analytical meanings 
is very hard. Hence we need to resort to numerical meth- 
ods. In particular, in this subsection we will numerically 
solve the evolution equations, derived using computer al- 
gebra, up to 10th order. An important question is then 
what values to choose for all the 63 moments involved at 
initial time <f>o = 0. Here we will simply choose the initial 
state corresponding to an unsqueezed Gaussian pulse in 
the volume, 

*(x) = ^ 7r ^e-^^ + — . (61) 

7T ' V ' 

While this state is very special in the space of all semi- 
classical moments, it already serves well to demonstrate 
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several features of quantum evolution. A more detailed 
analysis of the large parameter space will appear else- 
where. 

Using (|19j) . the initial moments are then given by, 

Min(o,&) 

ji r\ J — OO 

(62) 



d=0 



i — d 



where * denotes complex conjugation. This integral can 
be performed and happens to be vanishing except in 
those cases for which both a and b are even numbers, 



2—(a+b) foag.b—a a 



G a,b = 



! 61 



if a and b are even, 



(63) 







otherwise. 



This state saturat es the uncertainty relation, and we have 
that <t~ Ay = v 7 G 2 ' r~j 0(\fh), so the moments satisfy 
the semiclassical hierarchy. In our simulations, however, 
we have used % = 1 in order to make implications for the 
behavior of back-reaction more clearly visible. Neverthe- 
less, the features we point out have been verified also for 
smaller values of h. 

By this choice of initial values, the many-parameter 
family of all moments up to a given order is thus re- 
duced to a 1-parameter family labeled by the volume 
spread a. From this perspective, choosing a Gaussian 
initial state seems rather special and arbitrary; it is dis- 
tinguished only by the simplicity of writing down a wave 
function. As the system evolves, the state will depart 
from Gaussian behavior, a property which can easily be 
seen by following the behavior of the moments. Before 
discussing the numerical results in this direction, we note 
that the choice of Gaussian moments as initial values can 
sometimes be dangerous also for numerical purposes. A 
Gaussian state saturates the uncertainty relations, and 
numerical errors may easily move the evolved moments to 
the wrong side of saturation. Depending on the dynami- 
cal system, the saturation surface may be unstable such 
that uncertainty relations may be violated after some 
stretch of numerical evolution, which at that stage could 
no longer be trusted. (See also [TtJ for a discussion of 
similar issues in small- volume regimes.) In the numerical 
evolutions shown here, the fact that uncertainty relations 
are maintained has explicitly been tested by monitoring 
the relevant combinations of the second-order moments. 
In particular, the final time of each evolution has been 
chosen as the last time where the convergence tests are 
obeyed. 

Once an unsqueezed Gaussian is chosen as initial state, 
the only choice left is the value of a. The system under 
consideration has no ground state, and thus there is no 
preferred value for a as for instance the harmonic oscil- 
lator would suggest. For small initial fluctuations, which 
are related to a by a ~ AV ~ 1/AP, we just require 



that 1/Pq <C o <C Vq is satisfied for the initial values. 
Specifically, in all evolutions presented in this paper we 
take V = 1, P = 10 4 and a = 0.01. This does not 
represent a high restriction since we have also performed 
numerical simulations for different larger values of the 
width (cr = 1, 100, 10 4 ) and the qualitative picture is not 
changed. Even so, the behavior of the system for values 
of a < 10~ 4 is completely different. In this case the cor- 
rections to the classical trajectories are very large from 
the very beginning and the approximation can not be re- 
garded as valid. This happens because in the evolution 
equations ([25]) and flU, the moments G n ~ l x and G"'° 
appear. At the initial time </>o the former is zero, whereas 
the latter is of the order a~ n for even n and vanishing for 
odd n. Hence, for very small cr, the initial time deriva- 
tives of V and P increase much with the considered order. 
This entails a completely different trajectory for the clas- 
sical objects at each order and thus the approximation 
breaks down. 




4>div 
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FIG. 3: The ten different evolutions of the volume V at each 
order for A = 9 x 10 7 . The lowest trajectory corresponds to 
the classical one. The values are increasing with the consid- 
ered order in such a way that the largest trajectory is obtained 
at 10th order. For a more detailed analysis of the convergence, 
see App. [X] 



Finally, for the classical solutions (j3"5)) . (131)1) to exist, 
the cosmological constant must be in the interval [0, P§ ), 
so we have chosen three representative values: small A = 
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1, intermediate A = 10 4 , and large A = 9 x 10 7 . These 
values are much higher than the observed one. Even so, 
since the model we adopt as an example is unphysical 
anyway, it would not be meaningful to have it clad in 
a physical guise by using the observational value for A. 
We use the mentioned values in order to bring out more 
clearly the properties we have found about high-order 
moments. 



1. Numerical results 

Let us first explain the behavior of the volume V. For 
the case of small and intermediate values of the cosmolog- 
ical constant, the evolution of the volume reproduces the 
classical one with small corrections at all orders. How- 
ever, in the case of large cosmological constant we find 
that, at later times, the classical trajectory of the vol- 
ume receives large quantum corrections, which make it 
diverge faster (see Fig. [3]). In particular, we do not sec 
any indication that quantum back-reaction may trigger a 
recollapse; instead, the classical divergence of the volume 
is enhanced by quantum corrections with each order. In 
App. [X] we provide a complete analysis of the conver- 
gence of our solutions at different orders by studying the 
relative corrections that the expectation value of the vol- 
ume receives at each order. These results show that our 
effective analysis is valid (converges exponentially) at all 
considered times with a high precision. 

Regarding the evolution of the moments there are two 
different behaviors. On the one hand, those G lJ with 
(— l) l+ ^ = 1 only show one branch; see the two char- 
acteristic examples shown in Fig. 2) That is, their 
evolution at different orders reproduce the same trajec- 
tory (for large values of the cosmological constant with 
small corrections). On the other hand, the moments with 
(— iy+i = — 1 show two different branches. One of them 
always corresponds to the order i + j and the other one 
to higher orders. This curious behavior is due to the fact 
that, in the latter case, when truncating the equations 
of motion corresponding to this moment at order i + j, 
we neglect the contribution of certain moments that are 
initially nonvanishing. In the former case, for G l,J with 
(— 1) J+J = 1, all neglected moments in its equation of 
motion at order i + j happen to be initially vanishing. 

Before commenting on the collective behavior of dif- 
ferent moments, it will be useful to remember the par- 
ticular case with A = where, as shown in Sec. IIV Al 
analytic solutions for the moments at all orders can be 
obtained. In this case the evolution equations for the 
classical variables as well as for the quantum corrections 
are completely decoupled, and there is no quantum back- 
reaction. Using the exact solutions (|34p . the moments 
G ab with a > b are exponentially decreasing, whereas 
those with a < b increase exponentially. On the other 
hand, the variables with a = b are constants of mo- 
tion. Note also that the combinations \Z|G > l |/V and 
^/\G t '°\/P remain constant under evolution. In Fig. [5] we 




G 3 ^/n 3 
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FIG. 4: An example with two different moments with one 
and two branches respectively. For every moment G I,J , the 
evolution at every order from i+j to 10 is shown. The first two 
plots correspond to A = 10 4 , whereas the other two represent 
the evolution of the same moments for A = 9 x 10 7 . 
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FIG. 5: These figures show, in a logarithmic scale, objects 
that are constant in the case of a vanishing A. The plots 
correspond to the values A = 1, 10 4 , and 9 x 10 7 respectively. 
In the first two graphics the moments are constant throughout 
evolution until very late times. On the other hand, for a 
large cosmological constant this behavior is violated from the 
very beginning. As a general feature, in all cases, once t he 
moments are no longer constant, the quantities -\/\G '' \/V 
as well as \G l,l \ happen to be increasing functions, whereas 
\J\G l '°\/P are decreasing. The different lines correspond, from 
largest to smallest value at the initial time, to the moments 
G 4 ' 4 , G 2 ' 2 , and then, with an equal initial value by pairs, 
to (G 10 '°.G ' 10 ), (G 8 '°,G ' 8 ), (G 6 <°,G°> 6 ), (G 4 '°,G ' 4 ), and 
(G 2 '°,G ' 2 ). 



show these objects for different values of the cosmologi- 
cal constant. It can be clearly seen how the trajectories 
depart from the constant behavior earlier when increas- 
ing the value of the cosmological constant. In this simple 
system one can easily construct many other constants of 
motion, e.g. G a ' b / (P a V b ). Even though, for illustration, 
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FIG. 6: This figure shows the evolution of the moments G Z ' J 
with i > j and for both i and j even. Different plots corre- 
spond to the values A = 1, 10 4 , and 9 x I0 7 respectively. In 
the first two plots we observe that the exponentially decaying 
behavior of these moments is maintained for small and inter- 
mediate values of the cosmological constant. On the other 
hand, for large cosmological constant, the exponential decay 
is not followed anymore and there are even some moments 
that slightly increase. Interestingly, at the final time of the 
evolutions shown here, different moments G I,J with the same 
index i approach a common value. From largest to smallest 
values at the accumulation time the presented moments are: 
G 10 ' alone; G 8 ' 2 with G 8 ' ; G 6,4 with G 6 ' 2 and G 6 ' ; G 4 ' 2 
with G 4 ' ; and finally G 2 ' on its own. The tendency to this 
behavior can already be noted in the second plot, whereas in 
the first one all moments gather at a value near 1. 
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FIG. 7: This figure corresponds to G 1 ' 3 with i < j and both i 
and j even numbers. The value of the cosmological constant 
is again A = 1, 10 4 and 9 x 10 7 for different plots. As in the 
previous case, in the first two plots all moments coincide at 
a value near 1 at a given time, whereas in the last plot the 
different variables tend to gather in small sets given by those 
moments with the same number of P. In this last plot, at a 
final time and from largest to smallest value: G 4 ' 6 , then the 
three G 2 ' 1 , and finally the five G ' 1 . Spikes in this and the 
following plots arise from transitions of the moments through 
zero, shown in logarithmic form. 



log | G« 




log | G« 1/7^ 




A = 9 x 10 7 



FIG. 8: This figure corresponds to initially vanishing G IJ 
with i > j. The same values of A as in previous cases for 
each plot apply. Since the initial state we have picked is not 
adapted to the equations of motion, at the initial time all 
these modes are excited. Their absolute value after this ex- 
citation is approximately ordered by increasing (i — j). That 
is, the largest moment is G 9,0 , then G 9 ' 1 , G 8 ' 1 , ...On the 
other hand, it is also interesting to note that this excitation 
becomes larger with an increasing A. Finally, in the case with 
vanishing A, all these moments are exponentially decreasing. 
This behavior is qualitatively disappearing as one considers a 
larger cosmological constant. 



we will only consider those mentioned above. 

Finally, in Fig. [5] to UH we show the collective behavior 
of the different moments G 1 ' 3 . Essentially they have been 
divided in four groups depending on whether both i and 
j are even numbers or not and whether i < j or i > j. 
This classification is inherited on the one hand from our 
chosen initial conditions, where the only nonvanishing 
moments are those with both i and j even, and on the 



other hand from the case without cosmological constant 
explained above, where the increasing or decreasing be- 
havior of different moments also depends on the sign of 
i-j- 

The increasing behavior of those moments G 1,3 with 
i < j is obeyed also for nonvanishing values of the cos- 
mological constant (see Fig. [7] and [JJ]). On the other 
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FIG. 9: This figure shows the initially vanishing moments G l ' J 
with i < j. The same values of A as in previous cases for each 
plot apply. As in the previous case, initially all moments are 
excited to a given value, which is larger with an increasing 
(i — j). In this way, the moments with lowest absolute value 
are G ' 9 , G 1,9 , G 1 ' 8 , . . . On the other hand, and contrary to 
the previous case, the qualitative behavior of these moments 
is kept the same as in the case with A = even for large 
values of the cosmological constant. 



hand, in Fig. [5] and [5] the decreasing moments are plot- 
ted, but this behavior is violated as we move to higher 
values of the cosmological constant. As a general feature, 
in these plots it can also be observed that the evolution 
of different moments seems to be "slowed down" , that is, 
everything (for example the instant when nearly all mo- 



ments cross each other) happens latter as one increases 
the value of the cosmological constant. This makes the 
hierarchy between different moments be retained for a 
longer time. 

Another interesting issue is shown in the plots corre- 
sponding to the initially nonvanishing moments (Fig. [5] 
and [7]): at late times those moments G 1 ^ with the same 
number of P, that is with a common index i, accumulate 
at a common value. 

Finally, we note how the initially vanishing moments 
behave at the beginning of the simulations. As can be 
seen in Fig. [5] and |H1 since the initial state of the system 
is not adapted to the equations of motion, these moments 
are excited very quickly to a "natural" value which more 
closely corresponds to a dynamical coherent state. The 
absolute value of a given moment G 1,1 immediately 
after this excitation, increases with the value of the 
cosmological constant A and approximately also with the 
difference i — j. This last dependence can be seen in the 
equations of motion. More explicitly, making use of the 
general formula ([5]) and our particular Hamiltonian (|23[) 
it can be checked that the time derivative of a moment 
Qi,j i s given as a linear combination of the follow- 
ing objects {G i ' j - 1 G k - 1 '°,G l - 1 ' j G k - 10 ,G i ' j - 1 G k - 2 ' 1 

^Qi+k-nj-n^Qi+k-l-nJ+X-n}^ where fc mugt be 

summed from 2 to the order K at which we decide to 
truncate the Hamiltonian, and n corresponds to the sum 
over odd numbers that appears in Eq. ([5]). Taking into 
account the dependence of these objects on the Gaussian 
width a at the initial time, it is straightforward to 
see that the initial time derivative of G IJ behaves like 



2. Implications for state evolution 

In Sec. IIV CI several questions about the evolution of 
states has been suggested by the second-order analysis, 
but could not be answered without more general infor- 
mation about the high-order system of quantum back- 
reaction. With the numerical results, we can now provide 
additional indications regarding quantum evolution, but 
also new properties become visible. 

First, the strong increase of some of the moments al- 
ready seen at second order is confirmed. The system 
is thus necessarily one of strong quantum back-reaction. 
Results about self-adjoint extensions of quantum Hamil- 
tonians have suggested that the evolution can be contin- 
ued through the classical divergence of the volume. The 
only intuitive semiclassical interpretation would be that 
quantum corrections become so strong that they can trig- 
ger a recollapse. Our numerical results to the orders spec- 
ified do not provide any indication for this; they rather 
show that the divergence is enhanced by quantum cor- 
rections. 

The behavior of the moments confirms two expecta- 
tions about state evolution. First, several moments grow 
and become dominant providing strong quantum back- 



16 



reaction. Secondly, the state rapidly departs from the 
initially chosen Gaussian moments. While it is difficult 
to reconstruct a wave function from the moments, it is 
clear from several of their properties that the state can- 
not remain Gaussian. A Gaussian state has vanishing 
moments of odd order, a property which is immediately 
violated once the state is evolved. As the rapid initial 
increase of the odd-order moments shown in Figs. [5] and 
[9] shows, the evolution quickly adapts the state to one 
whose moments change less severely. To some degree, 
even the special choice of a Gaussian state may not be 
much of a restriction since the evolution soon leads to 
a more suitable dynamical state. However, for robust 
conclusions about quantum evolution from a fixed set of 
initial values one must analyze how the state settled down 
to depends on the initial moments. If different sets of ini- 
tial moments still lead to adapted states, but ones that 
differ from the Gaussian adapted state in a way sensitive 
to the initial set, the quantum behavior would be hard 
to predict without knowledge of what state may be pre- 
ferred. Here, much more numerical analysis of the large 
parameter space involved is necessary. 

Finally, we see from figures such as Figs. |5] and [9] that a 
hierarchy of the moments is maintained for a rather long 
time even throughout the phase of dynamically adapting 
the state. Moments of different orders clearly fall into 
distinct classes as for their behavior of decay or increase. 
Only when moments of different orders converge, which 
interestingly happens in a very narrow time frame for 
all orders shown, will the hierarchy be broken. At such 
a point, the state can no longer be considered semiclas- 
sical, and truncations of the infinite system of effective 
equations (and the asymptotic expansion they represent) 
become more difficult to justify in general terms. Nev- 
ertheless, in our specific case, even when the state is no 
longer of recognizably semiclassical form, the truncation 
is still justified by the convergence analysis we perform 
in App. [3] We have also shown the evolution even after 
that time of accumulation because it indicates another 
feature. After the accumulation, the moments again sep- 
arate into clearly demarcated sets, indicating that an- 
other hierarchy arises. However, since one has to evolve 
through an anhierarchical point, it is not clear whether 
this feature is one of the full system or an artifact of 
the truncation. Finally and somewhat unexpectedly, the 
hierarchy is maintained longer for larger values of the 
cosmological constant, even though the deviation from 
the harmonic model is then stronger. 



V. DISCUSSION 

The main contribution of this article is the introduc- 
tion of a new computational method to analyze quantum 
back-reaction of quantum mechanical or quantum cosmo- 
logical systems. We have shown that the use of efficient 
computer algebra tools in combination with the closed 
formula for the Poisson brackets of two generic moments 



has been essential to push the feasibility of computations 
to very high order. In particular, the example we have 
studied of a spatially flat, isotropic universe with a posi- 
tive cosmological constant and a free, massless scalar field 
already indicates the usefulness of these methods. 

Our analysis has found several new properties of state 
evolution, some of which were quite unexpected. For ex- 
ample, the state rapidly deviates from the initial Gaus- 
sian form (in the volume), but then settles down to an- 
other shape obeying a hierarchy of the moments. In this 
range, truncations used to analyze effective equations re- 
main justified. At some point, the moments converge, 
interestingly at about the same time. Several proper- 
ties found here remain without an analytical explana- 
tion, stimulating further studies. For instance, somewhat 
countcrintuitively, the moments happen to "slow down" 
their evolution as the value of the cosmological constant 
increases. That is, qualitatively they follow the same pat- 
tern, independently of A, but everything happens later 
for large values of the cosmological constant. Another 
interesting feature comes from the evolution of the ini- 
tially non- vanishing moments G 1 * 3 : at late times they 
accumulate at a common value for each index i. 

On the other hand, we have also studied the conver- 
gence of this truncated system of equations with respect 
to the order (see App. [SJ by analyzing the relative er- 
ror in the trajectory of the expectation value of the vol- 
ume V. Remarkably, even though a priori we expect 
only asymptotic rather than convergent expansions, the 
results converge exponentially for all considered times 
within a large range of orders. This means that our re- 
sults reliably reproduce the full quantum behavior of the 
system even quite near the divergence. 

Finally, details of the numerics remain to be explored, 
most importantly those related to the large parameter 
space involved. For instance, it is not clear yet how strong 
the role of the choice of an initial state is. The shape of 
the state changes rapidly in a very brief initial phase, 
as shown by a large change in the moments, and then 
settles down to a form better conserved by the dynam- 
ics. This evolved state seems adapted to the dynamics, 
but it is not known at present whether differently chosen 
initial states will give rise to the same kind of dynami- 
cally evolving state, nor is it known whether the initial 
choice could influence the dynamics strongly. For such 
questions, the parametrization by moments, rather than 
wave functions, is important because it gives full access 
to the state space. For instance, one could use a random 
number generator to construct the initial moments just 
by restricting to those sets that obey the Schwartz in- 
equalities. This will provide a systematic control to map 
the whole state space, which would not be achievable by 
specifying explicit wave functions. (Of course, one could 
randomize the coefficients of wave functions in some ba- 
sis, but the observable meaning of those variations would 
be much less clear than changes of moments.) Probing 
the large parameter space of the initial state is the main 
problem in this context which will benefit from further 
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numerical support. The results of this paper thus show 
that quantum cosmology provides its own set of prob- 
lems which are interesting from a numerical perspective, 
whose solution will then give feedback for the specific 
form of the dynamics realized. 
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Appendix A: Convergence of the system with an 
increasing number of moments 

In this appendix we address the issue of the conver- 
gence of the solution when considering an increasing 
number of moments. As we have explained in the main 
body of the article, the existence of a hierarchy on the 
moments define a semiclassical regime where the trun- 
cation of the system into a finite number of moments 
makes sense. But, since we have obtained the numerical 
solution of the system of equations at different orders, 
another method is at our hand to check whether the ig- 
nored moments are indeed negligible, namely to analyze 
if the difference between the expectation values at differ- 
ent order decreases sufficiently fast. 

In order to do so, we define the relative error at order 
n as, 



log AV n 
♦ 



1 



v n 



+1 



(Al) 



where V n is the expectation value at order n. If this 
object is convergent sufficiently fast with n, it gives an 
estimate of the total relative error (1 — V^/Kx,) commit- 
ted when truncating the system at a given order n. 

Even though we did not write it explicitly, the object 
AV n is time (<fi) dependent and, in fact, we expect it to 
increase (and eventually not to converge) as we approach 
the regime of large moments. Hence we have chosen three 
different times to perform the convergence tests: 0.4<^>fi n , 

0. 60fi n , and 0.9</>fi n , where </>fi n is the final time of each 
numerical evolution and takes the value 0.893</)di v for A = 

1, O.8O70div for A = 10 4 , and O.6650 d iv for A = 9 x 10 7 . 
These convergence tests are shown in Fig. [TJJ]in a nat- 
ural logarithmic scale. Firstly, we note that the conver- 
gence is exponential in all the cases and for all considered 
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FIG. 10: In these plots we show the relative change of the 
volume AV n at each perturbative order for the three differ- 
ent values of the cosmological constant in a logarithmic scale. 
The computed results and their linear regressions at times 
4> — O.40fl n , <j> = O.60fl n , and <j> — 0.9</>fi n are plotted respec- 
tively in black (dots and dot-dashed line), green (squares and 
continuous line) and red (diamonds and dashed line). In the 
first plot the points at n = 2 corresponding to <f> = O.40fl n and 
4> = 0.6</>fln are missing because they are exactly (up to our 
numerical error) zero. These two points have not been con- 
sidered to perform the corresponding linear regression. The 
slope of each linear regression gives the convergence order. 
The absolute value of the different slopes are, from upper to 
the lower plot and line: 1.22, 1.14, 1.02; 1.19, 0.98, 0.92; and 
0.59, 0.87, 0.99. 
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times. Therefore, we have performed linear regressions 
for all the data. 

As expected, the magnitude of the errors increases for 
later times in all cases. Even so, surprisingly for the first 
two plots, the convergence is faster at those times. The 
numerical value of the errors is also increasing with the 
value of A. For instance, the relative errors of the classical 
solutions (n = 1) are of the order e _6 — e -4 ~ 10~ 3 — 10~ 2 
for large A = 9 x 10 7 , whereas for small A = 1 they are 
only of order e~ 16 — e~ n s» 10~ 7 — 10~ 5 . This gives an 
idea of the magnitude of the back reaction in each case. 

On the other hand, we sec that our result at 10th order 
mimics very accurately the behavior of the whole quan- 
tum system, in the sense of an asymptotic expansion. In 
particular, the largest error we find corresponds to the 
large cosmological constant case at late times (O.90g n ) 
and it is of the order of e -9 w 1CP 4 . 

Note that in the last plot of Fig. [TO] (and also in the sec- 
ond one, but less severe so) even and odd orders show dif- 
ferent convergence behaviors. This observation demon- 
strates that moments of odd orders continue to contribute 
less significantly to the volume expectation value than 
moments of even order, even after the state has evolved 
away from Gaussian form for which odd-order moments 
vanish. [For odd n, the expectation values in the ratio 
V n /V n +i in (|A1[) differ by even-order moments, and the 
errors are seen to be enlarged in the plots.] 

This analysis proves that our treatment provides a 
valid approximation at all considered times, a result 



which strengthens the motivation to study this general 
approach. 



Appendix B: High-order equations 

In order to give a flavor of the increase in complexity 
of the equations at high orders, in Fig. [IT] we show the 
average number of terms per equation at each order. In 
addition, in this appendix, we also present explicitly the 
complete fifth-order equations of motion (with h = 1). 

Average number of terms 
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FIG. 11: In this plot we show the average number of terms 
that appear in each evolution equation with respect to the 
order. The dependence is exponential. 
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V 2 (P 2 - A) 5/2 8 (P 2 - A) 7/2 J 2 (P 2 - A) 5/2 

£ 2 3 _ ^,1 / 45PA 27PAG 2,2 \ + Q ? n f 3G 13 A 27PVG 2 ' 2 A 9PVA \ 



5/2 



8(P 2 -A) 5/2 4(P 2 -A) 5/2 / V 2(P 2 -A) 3/2 4(P 2 -A) 5/2 8 (P 2 - A) 

u / 9AG 2 ' 2 9A \ 3PG 2 - 3 9FA (4P 2 + A) G 3 - 9A (4P 2 + A) G 3 - 1 

V2(P 2 -A) 3/2 4(P 2 -A) 3/2 / 2VP 2 - A 8(P 2 -A) 7/2 4(P 2 -A) 7/2 
9FAG 3 - 2 45PVA (4P 2 + 3A) G 4 <° Ql i ( 3PAG 3 ' _ 3A (4P 2 + A) G 40 \ 



2 (P 2 - A) 3/2 16 (P 2 - A) 9/2 V 2 (P 2 - A) 5/2 8 (P 2 - A) 7/2 / 

G 2.2 ( 9VA ( 4p2 + A ) g3 '° 9A ( 4p2 + A ) ° 31 A5PVA (AP 2 + 3A) G 4 ' \ 
V 4(P 2 -A) 7/2 4(P 2 -A) 7/2 16(P 2 -A) 9/2 J 

9A(G 2 > ) 2 ( 9PAG 3,0 9A (4P 2 + A) G 4 <° \ 2 9PG 3,0 9AG 4 <° 9PAG 5,0 



£3,0 = ai H^ ; + ^ _ ^ V" ^ J v ^ G 2 : o _ ^ J ^ + 

4(P 2 -A) 3/2 V4(P 2 -A) 5/2 16(P 2 -A) 7/2 J 2VP 2 - A 4(P 2 -A) 3/2 4(P 2 -A) 5/2 

£ 3 ! _ WA (AP 2 + A) (G 3 ' ) 2 3AG 1 ' 1 G 3 ( 3A (4P 2 + A) G 31 15PVA (AP 2 + 3A) G 4 '° \ 3 
4(P 2 -A) 7/2 2(P 2 -A) 3/2 \ 4(P 2 -A) 7/2 16(P 2 -A) 9/2 / 



+ G 



, „ , 9AG 2 ' 1 9PyAG 3 '° \ _ 3PG 3 ' 1 3FAG 4 ' _ 9 A (4P 2 + A) G 21 G 4 >° 

\ 4(P 2 - A) 3/2 4(P 2 - A) 5/2 / VP 2 - A 2 (P 2 - A) 3/2 16 (P 2 - A) 7/2 

AG 4 ' 1 9PFAG 5 ' 



5/2 



4(P 2 -A) 3/2 4(P 2 -A) 

£ 3i2 _ 3A (4P 2 + A) (G 3 ' 1 ) 2 3AG 1 ' 1 G 3 ' 1 / 9PAG 2 ' 1 15PFA (4P 2 + 3A) G 4 \ q31 
2(P 2 -A) 7/2 (P 2 -A) 3/2 \ 2(P 2 -A) 5/2 8(P 2 -A) 9/2 / 

g2 /_ 9G 2 - 2 A 9PVG 31 A 9A \ ^ / W (AP 2 + A) G^A 27PA 



4(P 2 -A) 3/2 2(P 2 -A) 5/2 8(P 2 -A) 3/2 / V 2(P 2 -A) 7/2 8(P 2 -A) 5/2 / 

3PG 3 < 2 27A (4P 2 + A) G 4 '° 2 / 9PAG 3 ' 9A (4P 2 + A) G 4 -° \ 3FAG 4,1 



2VP 2 - A 16(P 2 -A) 7/2 ^ V4(P 2 -A) 5/2 16(P 2 -A) 7/2 / (P 2 - A) 3/2 
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3PA(G 3 ' ) 2 3AG 2 '°G 3 -° 3A (4P 2 + A) G 4 '°G 3 >° 6PG 4 * 3AG 5 -° 
(P 2 -A) 5/2 (P 2 -A) 3/2 4(P 2 -A) 7/2 VP 2 - A (p2_ A )3/2 

15PVA (4P 2 + 3A) (G 4,0 ) 2 3AG 1 > 1 G 4 -° QPAG^G 4 - 3^A (4P 2 + A) G™G A ZPAG^G 3 ' 1 

16(P 2 -A) 9/2 2(P 2 -A) 3/2 4(P 2 -A) 5/2 4(P 2 -A) 7/2 (P 2 - A) 5/2 

G 2 / 3AG 3 ' 1 9PyAG 4 '° \ _ 9PG 4,1 _ 3FAG 5 ' 

V (P 2 - A) 3/2 4 (P 2 - A) 5/2 ) ~ 2^/P 2- 1 _ 2 (P 2 - A) 3/2 

15A (4P 2 + A) (G 4 - ) 2 15AG 2 '°G 40 15PAG 3 -°G 4 '° 15PG 5 - 
16(P 2 -A) 7/2 4(P 2 -A) 3/2 4(P 2 -A) 5/2 2VP 2 - A 
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